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ABSTRACT 

A description of strange quark matter at finite temperature is given, with empha- 
sis on the inclusion of finite size effects. For massless quarks all thermodynamic 
quantities can be evaluated analytically, while the case of massive quarks must be 
treated numerically. The resulting model is applied to the phase equilibrium be- 
tween strange hadronic matter and strangelets, and some preliminary results are 
given. 

1. Introduction 

Low baryon number droplets of strange quark matter, so-called strangelets, 
might condense out of a hadronizing quark-gluon plasma in an ultra-relativistic 
heavy ion collision. The study of the ground state properties of strangelets is impor- 
tant in order to determine the possible (meta-) stability of these objects. However, 
even if strangelets are absolutely stable in the experimentally accessible region they 
may never be formed under the prevailing experimental conditions. Under which 
conditions we might expect to see strangelets cannot be inferred from a knowledge 
of the ground state properties of strangelets alone. In order to gain insight into 
the possible production of strangelets it is necessary to include the excited states 
through which the strangelet would pass at high temperature before ending up in its 
ground state. Using the MIT bag model^ for the strangelet the strong interaction 
is mimicked by the bag constant and the bag boundary conditions. Further inter- 
actions can be included in perturbation theory, but we have set = throughout. 
The single particle energy levels can be found by solving the Dirac equation with 
bag boundary conditions numerically. This is the shell model (see Madsen elsewhere 
in these proceeedings). An alternative approach is to compute a continous density 
of states based on the boundary conditions, which gives the liquid drop model. 
The latter approach is computationally much less expensive, and allows analytic 
evaluation of some quantities. Either of these models can be used to construct a 
partition function used for a statistical mechanical/thermodynamical description 
of strangelets. Shell model calculations have been carried out at finite tempera- 
ture by Mardor and Svetitsky^, but only in the case of vanishing baryon chemical 



potential. In the present paper only the liquid drop approach will be treated. A 
thorough study of the phase equilibrium between strange quark matter and strange 
hadronic matter was carried out by Lee and Heinz"^. They also calculated isentropic 
expansion trajectories for a number of different initial conditions. This work is an 
attempt to extend their treatment to include the important finite size effects in the 
quark-gluon phase. 

2. General Framework 

In this section we describe the general framework for a statistical treatment of 
non-interacting quarks and gluons in an MIT bag, using the liquid drop model. 

The starting point for the description is the grand canonical partition function, 
Z, for a collection of quarks, anti-quarks, and gluons obeying the confining boundary 
conditions of the MIT bag, which is considered to be spherical. In units where 
h = c = ks = 1^ Z IS given by: 

log^ = EE6-log(l - ^BV. (1) 

i k, 

B is the energy density of the perturbative vacuum inside the bag, known as the 
bag constant, and /3 is the inverse temperature. The index i runs over the species 
present in the bag. is -|-1 for fermions and —1 for bosons. The sum runs over all 
single particle states. In the liquid drop model this sum is replaced by an integral, 
by introducing the density of states, Pi{k), which as a consequence of the assumed 
spherical symmetry only depends on the length of the vector k. 

logZ = ^6- r dkp,{k)\og{l + (,e-P(V^-^^))- [iBV. (2) 

By applying a multiple reflection expansion'*'^ the density of states can be cast into 
the general form: 

"■('■) = {sj'-'^' + (f) + 1^ (t) ^' + ■ ■ ■} ■ '^') 

Here V = |7ri?"^ is the volume of the bag, S = AirR^ is the surface area, and C = SirR 
is the extrinsic curvature of the bag surface. The overall factor gi is the degeneracy 
related to the spin and colour degrees of freedom (6 for quarks and anti-quarks, and 
16 for gluons). The volume term is seen to be universal, i.e. it depends neither on 
the particle species nor on the imposed boundary conditions. The contrary is true 
for the surface and curvature terms. 



The surface term for quarks and anti-quarks has been known for some time^ to 




which vanishes in the hmit niq = 0. The surface term for gluons is zero. Thus for 
massless quarks, anti-quarks, and gluons the leading correction to the density of 
states is the curvature term. 

The curvature term for massive quarks and anti-quarks, has not been derived 
using the multiple reflection expansion method, but recently Madsen'' proposed the 
ansatz 

■'^ \ k J 127r2 \ 2mg \2 rnjy ^ ' 

which has the correct limiting values for = and — )• oo. For gluons the 
curvature term is® 

All thermodynamic quantities are obtained through the (grand canonical) ther- 
modynamic potential, 0, which, neglecting all higher order terms in Eq. (3) can be 
written 

= -riogZ 

= iii)y + 0,, 5(r, ii,)s + o,, a(r, ii,)C} + bv. (7) 

i 

In the case of zero temperature, the resulting expressions for can be found in 
the contribution by Madsen to these proceedings. 

For non-zero temperature we have not found an analytic expression valid in 
general. In the special case of zero quark mass, however, it is possible to evaluate the 
combinations Og^ + Ogj, for a single quark/anti-quark species, under the assumption 
/ig = — /Ug, which is valid when processes such as qq 77 are in equilibrium. 

^.c + n,,c = ^^^ + ^/^^ (9) 

The number densities are obtained as = — which for zero quark mass 

gives 

nq,v - riqy = i^gT^ + ^1^1 (10) 



1 , , 

ng,c-ng-^c = (llj 



The gluon contribution to the thermodynamic potential evaluates to: 

^.v = -^T^ (12) 
= \t\ (13) 



The free energy is obtained as 

F = n-J2N,iJ^,, (14) 

where Nq is the number of quarks minus the number of anti-quarks of flavour q: 

Nq = [nqy - nqy)V + [Uq^S " nq^s)S + [nqfi- - nqfi)C. (15) 

3. Quark Matter — Hadronic Matter Equilibrium 

In this section we describe how to incorporate flnite size effects into the treat- 
ment of phase equihbrium between quark matter and hadronic matter, at flnite 
temperature. 

3.1. Scenario 

In the context of heavy ion coUisions the formation of a strangelet wiU require 
a strangeness enriched plasma. If a quark-gluon plasma has a sufflcient strangeness 
content it may cool to form a (meta-) stable strangelet. 

Strange quarks are produced signiflcantly only in strong interaction processes, 
i.e. as part of ss-pairs. In order to develop a net strangeness the plasma must rid 
itself of some of the s-quarks during the process of hadronization. This is thought to 
happen through K^^K'^ emission^. We consider the resulting system, consisting of 
two phases, as having a flxed strangeness content, and as being in thermal, chemical, 
and mechanical equilibrium. This ideal situation is clearly not realized in heavy ion 
collisions, but serves as a reference for more soflsticated approaches as well as a 
testing ground for new ideas and models. 

3.2. The Model 

The quark-gluon phase is modelled as a collection of non-interacting u^d^s 
quarks, their anti-quarks, and gluons in a spherical MIT bag. 

A common u, d mass, m^, and chemical potential, /ig, is assumed. Collectively u 
and d will be referred to as the light quarks. There are thus only two independent 
chemical potentials, fiq and /i^, or fiB and /i5, where the latter are the chemical po- 
tentials for baryon number and strangeness respectively. The connections between 
these are 

11 
P'g = -^P'B P's = -P'B - P'S- (16) 

The hadronic phase consists of the known hadronic spectrum, as obtained from 
the Particle Data Group^° summary tables. Interactions are included by means 
of a hard core repulsion, leading to an excluded volume effect. The longer range 
attractive part of the potential has been neglected. There are several ways of 



including excluded volume efFects^^, but we use the statistical bootstrap method of 
Hagedorn and Rafelski^^. In this model thermodynamic quantities such as pressure, 
energy density, and baryon number density are all obtained from the corresponding 
ideal gas quantities by multiplication with the factor (1 + e(°)/(4i?))~^ , where the 
superscript (0) denotes an ideal gas quantity, and e and B are the energy density 
and the bag constant. Thus the pressure in the hadronic phase is given by 

The corresponding point particle expressions are those of an ideal gas, given by 



e 



„.". . ^,„!r „r).|,f .A.-,-^, ,20) 

where i runs over the hadronic states, gi is the degeneracy factor of the species, 
and bi is the baryon number, which is ±1 for baryons/anti-baryons and zero for 
mesons. The signs in the denominator of the integrands correspond to Fermi-Dirac 
and Bose- Einstein statistics. 

3.3. Equilibrium Criteria 

The system consisting of a strangelet and a bulk hadronic phase is considered 
to be in chemical and thermal equilibrium, so the temperature and the strangeness- 
and baryon- chemical potentials are the same in the two phases 

r(^) = rw, = ,r = ,^s''. (21) 

Minimizing the free energy leads to the mechanical equilibrium condition 

J2 + 5 + -^^^,c) +B + p(^) = 0. (22) 

This has the interpretation that the pressure exerted by the quarks and gluons 
in the bag is countered by the bag pressure, i?, the hadronic pressure, and the 
surface tension and curvature force. Note that whereas the surface tension is always 
positive, and thus acts against the pressure, the curvature force can act both ways; 
depending on the value of the chemical potential. In equilibrium the chemical 
potential for a hadronic species, z, is given by its net light quark and strange quark 
content Pq^i and Ps,i as 

= Vq^.fiq + Vs,zl^s- (23) 



The equilibrium conditions Eq. (21) and Eq. (22) can then be solved numerically, 
yielding a set of triplets, (T, /ig,//^), defining the phase diagram. 



4. Results 

The quark phase is parametrized entirely in terms of the bag constant, i?, and 
the quark masses. In almost all cases we have fixed the strange quark mass at 
nis = 150 MeV, while u and d have been assumed massless in all calculations. 
The bag constant has been varied between B^^'^ = 145 MeV, where u^d matter is 
marginally stable in the bulk, and B^^'^ = 195 MeV, where u, J, s matter is unstable 
relative to a gas of A's, at zero temperature. 

The hadronic phase is characterized in terms of the particle spectrum used, and 
the value of the bag constant. In most cases we have only included the ground state 
baryons A, E, S, and A(1232), while the mesons included are tt, K and rj. This 
gives practically identical results compared to the case in which a full spectrum of 
all resonances below ~ 2000 MeV was used, except for B^^'^ well in excess of 200 
MeV. All the results depicted here are for the limited spectrum. 

The strangeness contents of the system is parametrized through the strangeness 
fraction, /s, defined as the ratio of minus the strangeness density to the baryon 
density 

/. ^ (24) 

The results obtained for bulk systems are identical to those obtained by Lee and 
Heinz^. 



4.1. The Phase Diagram at T = 

At zero temperature the quark thermodynamic potential is given by the analytic 
expressions in Madsen's contribution to these proceedings. The integrals occurring 
in the expressions for the hadronic properties can also be evaluated analytically. 
There are no non-strange mesons, and no anti-baryons. Solving the equilibrium 
conditions Eq. (21) and Eq. (22) gives the phase diagram in the fiq-fis plane as 
depicted in Fig. 1 for a bulk quark matter phase. The solid line is the coexistence 
line. Hadronic matter is the preferred phase inside the region bounded by the 
coexistence line, while the system is in the quark matter phase outside this region. 
Baryons are only present at T = if the chemical potential, /ii, of the species is 
greater than the mass, m^. This is the case in the region outside the dashed line. 
Each segment of the line is marked with a symbol denoting the particle species 
appearing when crossing the line from the inside. When the chemical potential for 
K^,K~ becomes equal to the kaon mass, the kaons would undergo a Bose- Einstein 
condensation. Since this cannot easily be included in the model the region of the 
phase diagram to the left of the line marked K is not accessible within the present 
framework. Kaons will be present to the right of this line only at finite temperature. 
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Figure 1: (a) Phase diagram for bulk strange quark matter and hadronic matter at T = 0. 
(b) The strangeness fraction in the quark (full line) and hadronic (dotted line) phases as a 
function of the light quark chemical potential. B = (190 MeV)"*, = 150 MeV. See text 
for further details. 




so at zero temperature the region inside the dashed line will be devoid of hadrons. 
Hadrons may thus be present in the region between the full and dashed curves. 
Coexistence curves for lower values of B will enclose a smaller region of the phase 
diagram. For values of i?^/'* ^ 165 MeV the coexistence curve will lie totally within 
the dashed curve, since quark matter will then be stable relative to a hadron gas. 
At the point where the phase coexistence line crosses the line where nucleons first 
appear there is a visible "kink" due to the fact that nucleons start to contribute to 
the pressure. The dashed-dotted line represents systems with a total strangeness 
fraction of 0.1. In the hadronic phase such a system would consist entirely of A^'s 
and A's, for the parameter values chosen. 

Fig. 1 (b) depicts the strangeness fraction in the two phases along the coexistence 
curve, parametrized by the light quark chemical potential fig. For a sufficiently high 
strangeness fraction there is no hadronic phase - with this choice of parameters - 
at zero temperature. 

4.2. The Phase Diagram for T ^ 

At finite temperature the quark matter phase contains thermal gluons and anti- 
quarks, while the hadronic phase contains non-strange mesons and anti-baryons. A 
phase diagram similar to the one in Fig. 1 can be constructed for a fixed tempera- 
ture. The effect of increasing the temperature is that the hadronic phase occupies a 
smaller region of the phase diagram. This is because there is a large number of mass- 
less degrees of freedom in the quark matter phase, whereas the lightest constituent 
in the hadronic phase is the pion with a mass of ~ 140 MeV. This means that the 



pressure increases faster in the quark phase than in the hadronic phase, favouring 
the former over the latter. From this argument it is also clear that an increase in 
the bag constant will have the opposite effect, since B effectively contributes to the 
hadronic pressure. The quark matter phase benefits from an addition of strange 
quarks, as opposed to the hadronic phase where the presence of strange hadrons 
tends to increase the free energy, because of their larger masses. Thus increasing 
the strangeness fraction will favour the quark phase. This effect is clearly seen in 
Fig. 2, where the phase diagram for a fixed strangeness fraction has been projected 
onto the ub-T plane. The solid curve is for a strangeness fraction of 0.1, whereas the 
dashed curve is for = 0. The area between two like curves corresponds to a mixed 
phase. The quark-mixed phase boundary is the curve at high baryon density, and 
the other curve is the boundary between the mixed phase and the hadronic phase. 

The effect of finite size of the quark phase is seen in Fig. 3, where phase diagrams 
for a bulk quark phase and a quark phase with baryon number A = 20 have been 
superimposed. The most notable difference is that the hadronic phase grows on 
account of the mixed phase when going from a bulk to a finite quark phase. 

5. Discussion 
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Figure 2: Phase diagram for a strangelet with baryon number A = 20, and a bulk hadronic 
phase. The full curve is for a strangeness fraction, = 0.1, while the dashed curve is for 
/, = 0.0 B^/'^ = 180 MeV, = 0, = 150 MeV. 
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Figure 3: Phase diagram for a bulk hadronic phase in equilibrium with a bulk quark phase 
(dashed line), and a quark phase with baryon number A = 20 (full line). B^^'^ = 165 MeV, 
THg = 0, = 150 MeV, fs = 0.1. 

We have described how to incorporate finite size effects into calculations on the 
equilibrium between strange quark matter and strange hadronic matter. The results 
displayed in Fig. 3 show a clear difference from the bulk case, but the detailed effects 
of this difference as a function of baryon number and other parameters remain to 
be investigated. The goal is to use such calculations to learn whether strangelets 
may be produced in heavy ion collisions. 

The work described here is still in progress. The model presented can be im- 
proved upon in several aspects. As an example the partition function given in 
Eq. (1) counts all states, including those that are not colour singlets. A partition 
function describing only colour singlet states can be constructed using the group 
theoretical projection technique developed by Redlich and Turko^"^. 
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